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\Q Abstract 

We present a new representation of solutions of the Benjamin-Ono equation that 
are periodic in space and time. Up to an additive constant and a Galilean transfor- 
mation, each of these solutions is a previously known, multi-periodic solution; however, 
the new representation unifies the subset of such solutions with a fixed spatial period 
and a continuously varying temporal period into a single network of smooth manifolds 
connected together by an infinite hierarchy of bifurcations. Our representation explic- 
itly describes the evolution of the Fourier modes of the solution as well as the particle 
trajectories in a mcromorphic representation of these solutions; therefore, we have also 
solved the problem of finding periodic solutions of the ordinary differential equation 
^j! governing these particles, including a description of a bifurcation mechanism for adding 

or removing particles without destroying periodicity. We illustrate the types of bifurca- 
tion that occur with several examples, including degenerate bifurcations not predicted 
by linearization about traveling waves. 
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1 Introduction 

The Benjamin-Ono equation is a model water wave equation for the propagation of unidirec- 
tional, weakly nonlinear internal waves in a deep, stratified fluid [7] El [19]. It is a non- linear, 
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non-local dispersive equation that, after a suitable choice of spatial and temporal scales, 
may be written 



Our motivation for studying time-periodic solutions of this equation was inspired by the 
analysis of Plotnikov, Toland and Iooss |21[ [T3] using the Nash-Moser implicit function 
theorem to prove the existence of non-trivial time periodic solutions of the two-dimensional 
water wave over an irrotational, incompressible, inviscid fluid. We hope to learn more 
about these solutions through direct numerical simulation. As a first step, in collaboration 
with D. Ambrose, the author has developed a numerical continuation method jH |3] for the 
computation of time-periodic solutions of non-linear PDE and used it to compute families 
of time-periodic solutions of the Benjamin-Ono equation, which shares many of the features 
of the water wave such as non-locality, but is much less expensive to compute. 

Because we came to this problem from the perspective of developing numerical tools 
that can also be used to study the full water wave equation, we did not take advantage of 
the existence of solitons or complete integr ability in our numerical study of the Benjamin- 
Ono equation. The purpose of the current paper is to bridge this connection, i.e. to show 
how the form of the exact solutions we deduced from numerical simulations is related to 
previously known, multi-periodic solutions [221 UHl US] • Our representation is quite different, 
describing time-periodic solutions in terms of the trajectories of the Fourier modes, which 
are expressed in terms of N particles (3j(t) moving through the unit disk of the complex 
plane. Thus, one of the main results of this paper is to show the relationship between 
the meromorphic solutions described e.g. in [8] and these multi-periodic solutions. Our 
representation also simplifies the computation and visualization of multi-periodic solutions. 
Rather than solve a system of non-linear algebraic equations at each x to find u{x) as 
was done in [03], we represent u{x) through its Fourier coefficients by finding the zeros 
/3j of a polynomial whose coefficients involve only a finite number of non-zero temporal 
Fourier modes. We find that plotting the trajectories of the particles /3j(t) often gives more 
information about the solution than making movies of u{x, t) directly. 

A key difference in our setup of the problem is that we wish to fix the spatial period 
once and for all (using e.g. 2ir) and describe all families of time-periodic solutions in which 
the temporal period depends continuously on the parameters of the family. This framework 
may be perceived as awkward and overly restrictive by some readers, and we agree that the 
most natural "periodic" generalization of the 7V-soliton solutions |15[ [T7] of an integrable 
system such as Benjamin-Ono are the iV-phase quasi-periodic (or multi-periodic) solutions 
[T2"l l22l [TU1 [T6] . However, our goal in this paper is not to study the behavior of these 






2 



solutions in the long wave-length limit, but rather to understand how all these families of 
solutions are connected (continuously) together through a hierarchy of bifurcations. The 
additional restriction of exact periodicity makes the bifurcation problem harder for inte- 
grable problems, but easier for other problems that can only be studied numerically. To 
our knowledge, bifurcation between levels in the hierarchy of multi-periodic solutions of 
Benjamin-Ono has not previously been discussed. Indeed, with the exception of [10], previ- 
ous representations of these solutions are missing a key degree of freedom, the mean, which 
must vary in order for these solutions to connect with each other. The most interesting 
result of this paper is that counting dimensions of nullspaces of the linearized problem does 
not predict certain degenerate bifurcations that allow for immediate jumps across several 
levels of the infinite hierarchy of time-periodic solutions. As a consequence, in our numerical 
studies, we found bifurcations from traveling waves to the second level of the hierarchy, and 
interior bifurcations from these solutions to the third level of the hierarchy, but never saw 
bifurcations from traveling waves directly to the third level of the hierarchy (as we did not 
know where to look for them). This will be important to keep in mind in problems such as 
the water wave, where exact solutions are not expected to be found. 

We believe we have accounted for all time-periodic solutions of the Benjamin-Ono equa- 
tion with a fixed spatial period, but do not know how to prove this. Even for the simplest 
case of a traveling wave, it is surprisingly difficult to prove that the solitary and periodic 
waves found by Benjamin [7] are the only possibilities; see [5] and Appendix |Bj For the 
closely related KdV equation [TJ [181 123] , substitution of u(x, t) = uq(x — ct) into the equa- 
tion leads to an ordinary differential equation for uq(x) with periodic solutions involving 
Jacobi elliptic functions; see e.g. [20J. However, for Benjamin-Ono, the equation for the 
traveling wave shape is non-local due to the Hilbert transform. Nevertheless, Amick and 
Toland [S] have shown that any traveling wave solution of Benjamin-Ono can be extended 
to the upper half-plane to agree with the real part of a bounded, holomorphic function 
satisfying a complex ODE; thus, in spite of non-locality, we are able to obtain uniqueness 
by solving an initial value problem for uo(x). Interestingly, these traveling wave shapes are 
rational functions of e lx , which are simpler than the cnoidal solutions of KdV. This analysis 
of traveling waves via holomorphic extension to the upper half-plane is similar in spirit 
to the analysis of rapidly decreasing solutions of Benjamin-Ono via the inverse scattering 
transform H3] . Thus, it may be possible to prove that we have accounted for all peri- 
odic solutions of Benjamin-Ono by developing a spatially periodic version of the 1ST that is 
analogous to the study of Bloch eigenfunctions and Riemann surfaces for the periodic KdV 
equation [THJ [6], but this has not yet been carried out. 

This paper is organized as follows. In Section [2j we describe meromorphic solutions [8] 
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of the Benjamin-Ono equation, which are a class of solutions represented by a system of N 
particles evolving in the upper half-plane (or, in our representation, the unit disk) according 
to a completely integrable ODE. We also show the relationship between the elementary 
symmetric functions of these particles and the Fourier coefficients of the solution, which 
were observed in numerical experiments to have trajectories in the complex plane consisting 
of epicycles involving a finite number of circular orbits. In Section [3j we summarize the 
results of j3] in the form of a theorem (not proved in |3j) enumerating all bifurcations from 
traveling waves to the second level of the hierarchy of time-periodic solutions. The new idea 
that allows us to prove the theorem is to show that all the zeros of a certain polynomial 
he inside the unit circle using Rouche's theorem. In Section [4] we state a theorem that 
parametrizes solutions at level M of the hierarchy through an explicit description of the 
particle trajectories Pj(t). This theorem also describes the way in which different levels of 
the hierarchy are connected together through bifurcation. The proof of this theorem shows 
the relationship with previous studies of multi-periodic solutions. In Section [5] we give 
several examples of degenerate and non-degenerate bifurcations between various levels of 
the hierarchy. We also use these examples to illustrate some of the topological changes that 
occur in the particle trajectories along the paths of solutions between bifurcation states. 
Finally, in Appendix [Aj we give a direct proof that the particles /3j(t) in our formulas lie 
inside the unit disk in the complex plane; in Appendix [Bj we discuss uniqueness of traveling 
wave solutions. 



2 Meromorphic Solutions and Particle Trajectories 

In this section, we consider spatially periodic solutions of the Benjamin-Ono equation, 

Ut = Hu xx — uu x . (2) 



Here H is the Hilbert transform defined in (jl|, which has the symbol H(k) = — isgn(fc). It 
is well known [8] that meromorphic solutions of ^ of the form 

[ N 2k ] 

u(x, i) = 2Re|X: eiklx+kt - xm _ i ) ( 3 ) 

exist, where A; is a real wave number and the particles xi(t) evolve in the upper half of the 
complex plane according to the equation 

^ = V .„ 2 " + f .„ 2 " , (KKN). (4) 

dt e~ lk { x m-xi) _ i e -ik{xi-x m ) — i \ — — / \ > 

m=l m=l 
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This representation is useful for studying the dynamics of solitons of ^ over M, which may 
be obtained from ^ in the long wave-length limit k — > 0. However, over a fixed periodic 
domain IR/27rZ, we have found it more convenient to work with particles @i(t) evolving in 
the unit disk of the complex plane, 

01 = e~ iSl £A:={z: \z\ < 1}, (5) 

where a bar denotes complex conjugation. Up to an additive constant and a Galilean 
transformation, the solution u in Q3j) may then be written 

N 

u(x,t) = a + ^2u Mt) (x), (6) 
1=1 

where we have included the mean ao as an additional parameter of the solution and up{x) 
is defined via 

4|/?|{cos(s-fl)-|/?|} , _ ie , 

^ ) = ITF1^M)' {p = me ] - (7) 

From Q or direct substitution into pj), the /?j are readily shown to satisfy 

A = E 1 Z ,-i + E + <a - »o)A, (i<i<JV). (8) 

m=1 P/ Pm m=1 Pi Pm 

It is awkward to work with up(x) in physical space; however, in Fourier space, it takes the 
simple form 

r 2/3l fc l, £;<0, 
u/3,k = < 0, k = 0, (9) 

2/? fc , fc > 0. 

As a result, the Fourier coefficients Cfc(t) of u(x,t) in ([6| are simply power sums of the 
particle trajectories, 

{an k = 0, 

2[pt(i) + ••• + /?*(*)], fc>0, 

where = c_fc for k < 0. In [3], it was found numerically that although the 0i often execute 
very complicated periodic orbits, the elementary symmetric functions Uj defined via 

ao = l, Oj= E Ph'--Piv (j = l,...,N) (11) 

h<-<lj 

have orbits that are circles (or at worst, epicycles involving a finite number of non-zero 
temporal Fourier coefficients) in the complex plane. As a consequence, the spatial Fourier 
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coefficients Ck = 2tr(£ fc ), where k > 1 and £ is the companion matrix of the polynomial 
P(z) = Yl(z — Pj), also have trajectories that are epicycles, which we noticed immediately 
in our numerical simulations. All the solutions in this paper will be of the form 

N N 

]> - w)\ = E^^'^iW^ = e ~^)> (12) 

1=1 j=0 

where 

N 

P(z,\) = J2(-l) i *iWz N - j (13) 

3=0 

is a monic polynomial in z with coefficients &j that are Laurent polynomials in A, and such 
that for any A on the unit circle S 1 in the complex plane, the roots . . . ,/3jv of P(-, A) lie 
inside the unit disk. 

We may express the solution u in (j6j directly in terms of P as follows: 

N N ( oo "| 

u{x,t) = a + J2 u m(x) = «o + E 4Re \ E^(*) fee ' te f ( 14 ) 
i=i i=i lfc=i J 



AW J I P(*,A) J VA = e 



z = e 



Note that Hu = 4Re{— i [{zdzP) / P — N]}. The choice z = e~ lx (as opposed to e +lx ) 
follows from the decision in ^ to have Fourier coefficients with positive wave numbers carry 
powers of (3 rather than j3, while the choice A = e~ tujt leads to a natural sign convention 



when we interpret the exponents in the Laurent polynomials <X/(A) in (13) as measures of 
the direction and velocity of the traveling waves obtained in certain limits. It was shown in 
j3] that ( 14 ) is a solution of Q if there is a constant 7 G R such that 

T-Poo^oo + ^00 [P20 + wPoi + (a - 4N)Pio] 

+ Poo [P20 + ujPoi + (ao - 4iY)Pio] + 2Pi Pio = 0, 



where 



p jk = (zd z y(xd x ) k p( z ,x) 



z=e 1J " ' 
A=e" i " t 



(16) 



The goal of this paper is to find explicit formulas for the solutions P(z, A) of (15), show 
how they fit in with the previously known families of multi-periodic solutions described in 
j22[ HOl [T6] , and determine how these families are connected together through bifurcation. 

3 Paths connecting arbitrary traveling waves 

In [3], a classification of bifurcations from traveling waves was proposed after all time- 
periodic solutions of the linearization of (j2j about traveling waves were found in closed form. 
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A numerical continuation method was then developed to follow paths of non-trivial time- 
periodic solutions beyond the realm of validity of the linearization until another traveling 
wave was reached (or until the solution blows up as the bifurcation parameter approaches 
a critical value). Through extensive data fitting of the numerical solutions, the exact form 
of the solutions on this path was deduced. In this section, we give an alternative formula 
for these exact solutions that unifies the three cases described in [3] and makes it possible 
to show that the roots j3\ of the polynomial P(-,X) are inside the unit circle for X £ S 1 . 

An TV- hump traveling wave is uniquely determined by the mean, ao, a complex number 
(5 G A, and a positive integer, N: 

N 1-3I/3I 2 
u ao , NyP {x 1 t) = a + Y^u m {x), (50)=^-^, C = a -N i _^ 2 . (17) 

Here Pi ranges over all iVth roots of (3. This solution may also be written 

1 - 3|/?| 2 

u ao ,N,p(x, t) = u N ^(x - ct) + c, un,p = N - - + Nufs(Nx), (18) 

where un,@{x) = Nui^Nx) is the iV-hump stationary solution; hence, the traveling wave 
moves to the right if c > 0. We can solve for c and ao in terms of the period T > and 
a speed index i/£Z indicating how many increments of t£ the wave moves to the right in 
one period: 

nn 2ttz/ 1-3I/3I 2 . , 

cT=— , «o = c + iV T -^. (19) 

In order to bifurcate to a non-trivial time-periodic solution, the period T must be related 
to an eigenvalue 



(n)(N-n), Kn<N-l, 

~ ~ (20) 



(n + l-iV) [ n + l + iv(l-i=^) 



n > N 



of the linear operator [3] governing the evolution of solutions of the linearization about the 
A^-hump stationary solution: 

2Trm (nv + NZ, l<n<N, 

UN,nT=—, l<me{ (21) 
iV [(n + l)u + NZ, n>N. 

This requirement on the oscillation index m enforces the condition that the linearized solu- 
tion over the stationary solution return to a phase shift of itself to account for the fact that 
the traveling wave has moved during this time; see [3]. Here we have used the fact that if 
u(x,t) = un^{x) is a stationary solution and 

U(x, t) = u(x - ct, t) + c (22) 
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is a traveling wave, then the solutions v and V of the linearizations about u and U, re- 
spectively, satisfy V(x,t) = v(x — ct,t). The parameter (3 together with the four integers 
(N, is, n, m) enumerate the bifurcations from traveling waves, which comprise the first level 
of the hierarchy of time-periodic solutions of the Benjamin-Ono equation, to the second 
level of this hierarchy. We will see later that other bifurcations from traveling waves to 
higher levels of the hierarchy also exist, which is interesting as they are not predicted from 
counting dimensions of nullspaces in the linearization. 

After (numerically) mapping out which bifurcations (TV, u, n, m) and (N' , v', n' , m!) were 
connected by paths of non-trivial solutions, it was found that N, N', v and u' can be chosen 
independently as long as 

N' 

N' < N, u'> —v. (23) 

The other parameters are then given by 

m = ml = Nv' — N'v > 0, n = N - N' , n' = N-l. (24) 

The following theorem proves that these numerical conjectures are correct. 

Theorem 1 Let N , N' , v and v 1 be integers satisfying N > N' > and m = Nv'-N'v > 0. 
There is a four-parameter family of time-periodic solutions connecting the traveling wave 
bifurcations (N 1 , u', N — l,m) and (N, v, N — N' , m). These solutions are of the form 



N 

u{x,t) = a + ^2up l{t _ to) (x - x ), (25) 
l=i 



where 0i(t), . . . , /?Ar(t) are the roots of the polynomial P(-, e tujt ) defined by 

N , aw'N-N' , n\v-v'N' 



P(z, A) = z N + A\ v z N ~ N + B\ v ~ u z N + C\ u , (26) 



with 



N . N i + s + s i (N + s')s' 

V N + s + s' V N'(N-N') + (N + s')s n 1 ' 



B _ I (N + s>)s' I s c _ f~ s /N-N' + s + s' 



N'(N-N') + (N + s')s']j N-N' + s' yN-N' + s\ N + s + s' ' 

N 2 v'-(N') 2 v 2N'(i/-v), 2vr N'(N — N')(N + 2s') 

a = 2s s, lo = — = . 

m m 1 m 

The four parameters are s > 0, s' > 0, xq G M and to £ R. The N- and N 1 -hump traveling 
waves occur when s' = and s = 0, respectively. When both are zero, we obtain the constant 
solution u(x,t) = — - - . 
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Proof: Without loss of generality, we may assume xq = and to = 0. It was shown in 



[3] that P(z,\) of the form (26) satisfies (15) if 



(28) 
(29) 



7 = (3iV - a )N - vuj, 

[(N' f - 2NN' + N'a - u'u)B + [(N') 2 + 2NN' - N'a + v'uAAC = 0, 
[3N 2 - ANN' + (N') 2 - (N - N')a + (u - BC 

- [N 2 - (N 1 ) 2 — (N — N')a + (v- v')u]A = 0, 
(Nolq -vu- N 2 ) + [(2N' - N)a + (y - 2i/)u + 3N 2 - 8NN' + A(N') 2 ] B 2 

+ [(N - 2N')a + A(N') 2 - N 2 + (2v - v)u] A 2 + [(3iV - a )N + uw] C 2 = 0. (31) 



(30) 



Using a computer algebra system, it is easy to check that (29)-(31) hold when A, B, C, oq 



and u are defined as in (27). When s' = 0, we have A = B = and C 

ojv N'(N - N')v 



-Ce 



-ict 



N 



m 



a - N 



N + s 

1-3C 2 
1-C 2 ' 



so that 



where each (3i is assigned a distinct iVth root of —C. By (17), this is an TV-hump traveling 
wave with speed index v and period T = Similarly, when s = 0, we have B = C = 

and A = 



— so that 

N' + s' 

Xf^Ae- ict I < N' 
I > N' 



uju' (N - N')(N + 2s'y .1-3.4 2 
c = ^rrr = = ao — N 



N' 



m 



I- A 2 ' 



which is an TV'-hump traveling wave with speed index v' and period T = 

Finally, we show that the roots of P(-, A) are inside the unit disk for any A on the unit 
circle, S . We will use Rouche's theorem [2]. Let 

= z k + a\»'z n ~ n ' + B\ v ~ v ' z N> + C\\ 



z N +A\ v 'z 

z N 



u'N-N' 



+ BX^'z N '. 



From ([27]), we see that {A, B, C} C [0, 1), ,4 > 5C, B > CA and C > ,45. Thus, 
:= \f2(z)\ 2 - \h(z) - f 2 (z)\ 2 = \\-»'z N ' + A\ 2 - \BX~ V ' z N ' + C\ 2 



1 + A 2 - B 2 - C 2 + 2(A - BC) cos 9 > (1 - A) 2 - (B - Cf 



(32) 



where A v z 



e l6 '. Similarly, 



^3 := \h{z)\ 2 - \h{z) - h(z)\ 2 > (1 - B) 2 -{A- C) 2 . 



(33) 
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Note that 



B < A, 


C <B = 


> B- 


C <B 


- AB < 1 


-A = 




> for z 6 S 1 


B < A, 


C > B = 


> \C- 


-A\<1 


- B 




> d 3 (z) 


> for z G S 1 


A< B, 


C <A = 


> A — 


C <A- 


- AB < 1 


- B = 


> M z ) 


> for z 6 S 1 


A< B, 


C> A = 


> \C- 


- B\ < 1 


-A 




> d 2 {z) 


> for z G S 1 



Thus, in all cases, fi(z) = P(z, A) has the same number of zeros inside S 1 as f2{z) or fs{z), 
which each have N roots inside S . Since is a polynomial of degree N, all the roots 

are inside S . □ 



4 An infinite hierarchy of interior bifurcations 

Next we wish to find all possible cascades of interior bifurcations from these already non- 
trivial solutions to more and more complicated time-periodic solutions. The most interesting 
consequence of the following theorem is that there are some traveling waves with more 
bifurcations to non-trivial time-periodic solutions than predicted by counting the dimension 
of the kernel of the linearization of the map measuring deviation from time-periodicity. This 
will be illustrated in various examples in Section [5] 

Theorem 2 Let M > 2 be an integer and choose 

ki, ■ ■ ■ , Um G N, (positive integers, not necessarily distinct or monotonic), 

k _i 

u\, . . . , vm £ (arbitrary integers satisfying Vj-\ < — — Uj for j > 2). 

kj 

Now define the positive quantities 

'", h ^ h", ■ rj = *i(*J + V-Ofcj-i , ^ (2<i<M). (34) 

Let J = argmax 2 <j<Af Tj. If there is a tie, J can be any of the candidates. Then there 
is an M + 2 parameter family of time-periodic solutions of the Benjamin-Ono equation 
parametrized by 

si > 0, sj > 0, x j0 G M, (1 < j < M) (35) 
and constructed as follows. First, we define 

Zl^Jl + ll Sj: (2<j<M, j^J), (36) 

Ql = si, pi = s 1 +k 1 , q j =p j _ 1 +s j , Pj =q j +k j , (2<j<M) (37) 
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so that Sj > for 1 < j < M and 

Si k\ S2 &2 S M-1 fcM-1 Sjvf fcjvf 

'^9i' X ^Pi / <^72 / <^' • •'^^/-r^^Af-l^^gA/'^PA/ • (38) 
For any subset S of Ai = {1, . . . , M}, we denote the complement by S' = M \ S and define 

ks = ^kj, u S ' = ^2 u m , Cs = I \\ a jm \i J] 6 m , (39) 
jes mes' V(i,m)e5xS' / VmeS' / 



where 

1 (Pm ~ qj)(q m ~ Pj) u _ [%n_ e -ik m x m0 _ / 4 qn 



Qj){Pm-Pj) 




Then 

N 



i(x,t) = a + J2 u Mt)( x ) ( 41 ) 



z=i 



is a periodic solution of the Benjamin- Ono equation, where N = kjy[ = X^jli ^ 

«o = ^2iV - A* + ^tj^ - 2si + ^Ijsj, w = ^ = r j + 7j sj, (42) 
and /3i(i), . . . , /3/v(i) are t/ie roots o/ i/ie polynomial z i— ► P(z, e~* wt ) given by 

P(z,X)= Y, C s \ Vs 'z k s. (43) 

When M = 2, i/iis representation coincides with that of Theorem^ if we set 

ki = N — N' , v\ = v — v' , s\ = s, xio = xo — -^uto, 

h 

fc 2 = N , v 2 = u', s 2 = s', x 2 o = x - —uit . 

k 2 

It reduces to a traveling wave when s = or s' = and to a constant solution when both 
are zero. Similarly, when M > 3, the solution reduces to a simpler solution in this same 
hierarchy (with M replaced by M = M — 1 ) when s\ or sj reaches zero. Specifically, if 
s\ = 0, then P(z, A) = z kl P(z, A), where P(z, A) corresponds to the parameters 

kj = kj+i, i>j = u j+1 , sj = Sj +1 , x j0 = Xj+i )Q , (1 < j < M). (44) 

We interpret this as an annihilation of k\ particles fl\ at the origin. If sj = 0, we have 
P(z, A) = P(z, A), where 



kj = kj, Vj = Uj, Sj = Sj, Xjo = Xjo, (1 < j < J — 2) 

kj = kj + kj, Vj = Uj + vj, Sj = Sj , x j0 = kjX i°X k k J J XJ ° , U = J~ 1)> 
kj = kj+i, Vj = v j+ i, Sj = s j+ i, Xj = Xj + i fl , (J <j< M). 
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If several Sj are zero when sj = (i.e. if a tie occurs when choosing J = argmax 2 <j<j^ Tj), 
the bifurcation is degenerate and any subset of the indices for which Sj = may be removed 
using the rule (45) repeatedly (once for each index removed, with J ranging over these 
indices in reverse order to avoid re-labeling), allowing for bifurcations that jump several 
levels in the hierarchy at once. 



Proof: Rather than show directly that P(z, A) in (43) satisfies (15), we show that each 
of our solutions differs from a multi-periodic solution [22, 16 j by at most a transformation 
of the form (22 ). We give a direct proof that all the zeros of P(-, A) lie inside the unit circle 
in Appendix |A| as concluding this from the combined results of ^U\ and |16] is complicated. 



(46) 



If all the inequalities in (38) are strict, then it is known [16J that 



U 



2i— log- 

OX J 



satisfies the Benjamin-Ono equation ^ with 




(M) 
j<m 



]m 



(M) 

+ ^ Hjfj, m Aj r 

j<m 



Pj-Qj, r i Pi ■ <U 



exp(A,- 



Qj){Pm Pj) ( c m c j) 2 (km kj 



(47) 

(48) 
(49) 
(50) 
indicates a 



(Pm - qj)(q m - Pj) (Cm - Cj) 2 - (k m + kj) 2 ' 

where Xjo £ M are M arbitrary phase parameters and the notation Y2u=o l 
summation over all possible combinations of /xi = 0,1; ^2 = 0,1; ...; = 0,1. (The 
notation Y^j< m indicates that j and m both vary between 1 and M such that j < m, while 



X^m^j indicates that m varies from 1 to M omitting m = j). 

We write /' and / in (47) and (48) as sums over all subsets S of A4 = {1, . . . , M}: 



/',/= E 

SeV(M) x (j,m)eSxS' 



n 



— lA- 

g 2 "jm 



n- 

its 



} ikj (x—Cjt—Xjo)^f(f>j 



(S' = M\S), 



(51) 



where —<j>j is used for f and +4>j is used for /. Next we observe that 



u(x,t) = U(x - ct,t) + c (52) 
will be time-periodic with period T = — if there exist integers Vj G Z such that 



kj (yC ~\~ Cj 



(1 < j <M). 



(53) 
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Then we have 



f( x -cM)=(n etkjx ) e c 4u e- i " mujt ) ( n e ~ ikjx ) > 

f(x-ct,t) = 



(54) 



SeV(M) 



•meS' 



n 



e ?5 n< 

5e"P(A4) V meS 



J65 



with C5 as in (39) above. The complex conjugation in / comes from interchanging S and 
5' in the sum after factoring out ( IljeX " ' ' )• ^ follows that u in (52) is given by 

'id x g 



u = c-2N + 2id x log f = a + 2 



AM +2 



-id r h 



N 



(55) 



where 



(56) 



AT = J>i, a = c + 2Af, 5 (x,t) = P(e- ix ,e- i " t ), h = g, 
j'=i 

S£V(M) m<=S' j€S 

With A fixed, P is a monic polynomial in z of degree A r . If we complexify x and fix t, then 



P(e 



— „— iu>t\ 



/'(x-ct,t) = 4^ f(x-ct,t)=Q, 



(57) 



so all the zeros of P are inside the unit circle iff all the zeros of / are in the upper half-plane 
and all the zeros of /' are in the lower half-plane. These properties of / and /' were assumed 
to be true in [22 , leaving a small gap in their proof (acknowledged in the paper); we give 
a proof in Appendix |Aj The right hand side of ( 55 ) is equal to the right hand side of ( 14 ) , 
which establishes the representation (41) of u(x,t) in terms of the trajectories @i(t) of the 
roots of P(-,e~ iujt ). 

Eliminating c from (53) and using (cj — Cj_i) = (Ay_i + kj + 2sj), we find that 

kjkj-i(kj-i + kj + 2s j) = (kj-ifj — kjUj-i)u, (2 < j < M). (58) 



This shows that rrij in (34) must be positive. Eliminating u>, we find that 
T j + = tj + Jjsj, j,J e {2, . . . , M}. 



(59) 



Choosing J = argmax 2 <,<j^ Tj and solving (59) for Sj in terms of sj yields (36), which 
ensures that Sj > whenever sj > 0. From 



a = c + 2N, 



-ci + y^u, ci = fci + 2si, uj = Tj + -fjsj, (60) 
fei 



we obtain the formulas in (42) for ao and ui. 
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Finally, we drop the assumption that the inequalities in ( 38 ) are strict and observe what 



happens to these solutions when si = or sj = 0. If si = 0, then b\ = 0, so 
1 €' 5 C s = 0, 1 G 5 ^ C s = ]"[ a jm )[ ]"[ 

\(i,m)e(S\{l})xS" / VmeS" 

But since q\ = and pi = /ci when s\ = 0, we have 



aimbm = (Prn-qi)(q m -pi) Jq^ e - lkmXm0 = Qrn^h e . lkmXm0 ^ ^ 
V {Qm ~ qi){Pm ~ Pi) V Pm V Pm ~ fa 



Thus, if we define M = M — 1 and shift indices down as in (44), the parameters qj and pj 
will decrease by fa as illustrated here, 

S1=0 ki S2 k-2 S3 fcs SM k M 

^=^Qi ^<^Pi ^<^q2 ^<^Pi < Q3 ^T^- • -^r^lM ^K^PM (62) 
O^J^Ql < Pi < <h < ^-\^ < -\ < PM-1 

S1=S 2 k!=k 2 s 2 =s 3 k 2= k 3 S M -1=SM k M _ 1= k M 

and hence dj m = aj + i )Tn+ i and b m = ai^ m+ ib m+ i for 1 < j, m < M, j ^ m. Therefore, with 
the notation 5' = M \ S, the only difference between 

P(z,X)= Yl CsX^'z^, P(z,X)= £ C^s'zh, 

seV(M),ies ser(M) 

is that each term in the former sum carries an extra factor of z kl when the sets S and S 
are matched up in the natural way. Thus P(z,A) = z kl P(z,\). 

Now consider the case sj = with J = argmax 2 <y <A/ / Tj. This time aj-ij = due to 
qj = pj-i, so Cs = unless J — 1 and J are both in S or both in S' . Let us define 

Vj(M) = {S £ V{M) : J-l,JeS or J-l,JeS'}. (63) 

When we perform the sum over S G T'j{M) to construct -P(z, A), we can consider J — 1 and 
J as a single unit. The following factors always appear together in any Cs that contains 
one of them: 



(Pm ~ qj-i)(q m ~ Pj) , , qj~i - 
aj-i, m aj m — \ t T7 r, — \ e 

V - 9J-l)(Pm - PJ) V W 

Thus we can remove pj-i and <?j from the sequence if we short circuit the diagram 

si fci S2 sj-i fcj-l sj=0 kj sj+i k J+1 s M k M 

< Qi < Pi < ■ • • < QJ-i < PJ-i = QJ < PJ < QJ+i < ■■■ < QM < PM 
< qi < PlJ^' < __^^ J_1 - ^ J < ' \3_J A/_1 < ^ M_1 

Sl fc i S2 SJ ~ 1 fcj_ 1 =fc J _ 1 +A:j sj=sj+l fej SM-1 k M -i_ 
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and set 

kj-i = fcj_i + kj, 



vj-i = vj-i + vj, 



kj-ix,j-i fi + kjxj 
Xj - lfi = Ft i+kj • (64) 



The other parameters are simply copied from the original sequence as was indicated in ( 45 ) . 
We then have 



P(z,X)= Yl C s \ u s'z k *= y, C ~s Xh ' 



z s 



P(z,X), 



sev.j(M) 



where again the sets S and S in these sums are in natural 1-1 correspondence. 

Finally, we verify that the parameters kj , Sj and Vj of the reduced system are consistent 
with the construction, i.e. if M > 2 and we define J = argmax 2< ^. < ^fj, then 



Tj-Tj 7? 



(2<j<M, j + J). 



Recall that these equations are obtained by eliminating c and u) from 

kj(c + Cj) = ujoj, (1 < j < M). 



In the first case where si = 0, (44) and (53) together with 

£ j = Pj + = c i+i ~ 2fe i> (1 < J < M) 



imply that ( 66 ) is satisfied if we define 

c = c + 2k\, 



uj = u. 



Since we annihilate k\ particles in this case, N = N — k\, and the mean 

a = c + 2N = c + 2N = a 



(65) 



(66) 



(67) 



(68) 



(69) 



remains unchanged in spite of the change in c and iV (as it must for the solution to vary 
continuously through the bifurcation). In the remaining case where sj = 0, we have 



Cj, j<J-l, 
Cj + kj = cj-kj, j = J-l, 
{Cj+i, J<j<M. 



(70) 



Thus, by (|45|) and (|53|), all the equations in (|66|), except possibly j = J — 1, are trivially 

(71) 



satisfied if we define 



c = c, 



UJ = Ul. 
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The remaining equation is 



(kj-% + kj)(c + cj_i + kj) = {vj-i + vj)uj. 



(72) 



To see that this is true, note that by (53), 



/cj_i(c + c/_i) = Vj-\U), 



kj(c + cj) 



UJUJ. 



(73) 



Adding these equations and using kjcj = kj{cj-\ + kj_i + kj) gives (72), as required. Since 
c = c and N = N, the mean ao = c + 2N does not change as a result of the bifurcation. 
Thus, we have shown that when si = or sj = 0, the function u(x,t) in (41) agrees 



with another function in the hierarchy with M reduced by one and the parameter si or sj 
removed. Continuing in this fashion, we can remove all the zero indices, eventually yielding 
the case where the pj and qj are distinct from one another, or the case that M € {0, 1}. This 
shows that u(x,t) is a solution of Q, where we rely on Theorem [T] to handle the reduction 
from M = 2 to the traveling wave case M = 1, or the constant solution case M = 0. □ 



5 Examples 

In this section we present several examples to illustrate the types of bifurcation that occur 
in the hierarchy of time-periodic solutions described in Theorem [2} We begin with the 
simplest example that leads to a degenerate bifurcation, namely 

M = 3, £=(1,1,1), v = (-2,-1,0). (74) 

This solution and the three M = 2 solutions connected to it have parameters shown in 
Figure [T] We hold the mean ao = 0.544375 fixed, which is the value used in several of the 
numerical simulations in |Jl |3] , and construct a single bifurcation diagram showing all four 
solutions; see Figure [2} Note that paths B,C,D are actually parametrized by 

-(B) -(B) -(C) -(C) -(D) -(D) ,- K s 

S\ = S 2 , 4 = S 3, S\ = Si, 4 = S 3, s\ = Si, s y 2 = S 2 (75) 



in formulas (34)-(43), but we use the original variables Sj here for the convenience of 
making a single bifurcation diagram. The one confusing aspect of doing this is that we 
obtain different traveling waves depending on the order in which we set the Sj to zero. This 
is why we drew two axes for s 2 and S3 in Figure [2] For example, if we start on path A and 
decrease s\ to zero (moving to path B) and then decrease S3 to zero, we obtain the traveling 
wave bifurcation (2, —1, 1, 1); however, if we first set S3 = (moving to path D) and then 
set si = 0, we obtain the bifurcation (2, —1, 2, 3). Both traveling waves have N = 2 humps 
and speed index v = — 1, but the amplitude and period of the two solutions are different as 
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u 
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v= -3,0, 
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= -2,-1, 


m = 


1,1, 


m 


1, 


m = 3, 


m 


3, 


T = 


2,2, 


T 


2, 


T= 2, 
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2,2, 
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2, 
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7 


= 4/3, 


S3 = 


■S2, 


Sl 


o, 


S2 = 0, 


S3 


o, 
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l-2si-4s 2 , 


Of) 


= l-2s 2 -2s 3 , 


a = l-2si-2s 3 , 
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= 1 -2si - -s 2 , 




2 + 2s 2 , 




= 2 + 2s 3 , 


4 




4 




uJ 


w =2 + -s 3 , 


uJ 


= 2+ 3* 



Figure 1: Parameters of four paths of time-periodic solutions connected by bifurcations. 

they have different bifurcation indices. Similarly, although starting on path A and setting 
■s 2 and s 3 to zero in either order leads to the same stationary solution, the period T is 
different depending on whether we follow path B to (1,0, 1, 1) or path C to (1,0, 2, 3). This 
only happens at the bottom (M = 1) level of the hierarchy, where k%, v\ and s\ are not 
sufficient to uniquely determine the traveling wave; if we start with M > 4 and follow two 
paths down several levels to M > 2 with the same parameters kj, Vj and Sj, the resulting 
solution is independent of the path. 

In Figure [3] we plot the particle trajectories of several solutions on paths C, D and A 
in the bifurcation diagram of Figure [2] We parametrize each path linearly by a variable 
9 £ [0,1]. For example, on path A, 

Sl = l -^e, s 2 = sz= l -^{i-e), o<0<i. (76) 

Path C connects the one-hump stationary solution to the three-hump traveling wave. When 
9 = 4.0 x 10 -6 , two particles have nucleated at the origin and execute small, nearly circular 
orbits around each other while the third particle travels around its original resting position. 
As 9 increases, the orbits deform and coalesce into a single path, as shown in the middle 
two panels of this row. At the critical value 9 = 0.002649485, the particles collide at 
t = T/6, t = 3T/6 and t = 5T/6, so the solution of the ODE ([8| ceases to exist for all 
time; nevertheless, u(x,t) in (41) remains smooth and satisfies Q for all t. As 9 increases 



to 1, the common trajectory of the three particles becomes more and more circular until 
the traveling wave is reached, where it is exactly circular. The solutions on this path are 
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Figure 2: Left: degenerate bifurcation from M = 1 to M = 2 (paths C and D) and M = 3 
(path A) with ckq held fixed. Right: three dimensional plot of the solution labeled D4. 
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Figure 3: Particle trajectories along paths C,D,A in the bifurcation diagram of Figure [2] 
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reducible in the sense that their natural period is 1 /3 of the period T used here. However, 
in order to bifurcate to paths A and D, we have to use this solution rather than the reduced 
solution. 

Path D connects the two-hump traveling wave with speed index v = — 1 to the three- 
hump traveling wave with speed index v = —3. When we bifurcate from the two-hump 
traveling wave, a new particle nucleates at the origin and grows in amplitude until its 
trajectory joins up with the orbits of the outer particles. As increases further, the three 
orbits become nearly circular and eventually coalesce into a single circular orbit at the 
three-hump traveling wave. Note that the particles on path D follow different trajectories, 
which leads to a "braided" effect in the peaks and troughs of the solution u(x, t) shown in 
Figure [2j unlike path C, these solutions are not reducible to a shorter period. 

Path A connects an interior bifurcation on path B to this same three-hump traveling 
wave. When 9 = 0.0005, a particle has nucleated at the origin without destroying the 
periodicity of the orbit of the other two particles. This path involves two topological changes 
in the particle trajectories, as shown in the middle two panels of the bottom row of Figure|3| 
As 6 — > 1, this solution also approaches the three- hump traveling wave. This is interesting 
because, up to a phase shift in space and time, the linearized Benjamin-Ono equation 
about this traveling wave [3] has only two linearly independent, time-periodic solutions 
corresponding to the bifurcations (3,-3,1,3) and (3,-3,2,3); this degenerate bifurcation 
is not predicted by linear theory. 

In a similar way, we can construct a bifurcation from a traveling wave to an arbitrary 
level of the hierarchy by taking 

M arbitrary, £=(1,1,..., 1), v = (— M, —M + 1, . . . , —2, —1, 0). (77) 

We find that m,- = 1, Tj = 2 and jj = 2 for 2 < j < M, so any subset of the indices 
J = 2, . . . , M can be removed to obtain a solution at a lower level of the hierarchy. If all 
the indices are removed, we obtain a traveling wave with 2 M_1 — 1 bifurcations to higher 
levels of the hierarchy, but only M — 1 of them (to the second level) are predicted by linear 
theory. The case M = 4 is depicted in the bifurcation diagram of Figure |4j where each 
tetrahedron contains paths with one of the Sj set to zero, and the outermost point of the 
outer three tetrahedra corresponds to one and the same traveling wave. Linear theory 
predicts the bifurcations (4,-6,1,6), (4,-6,2,8), (4,-6,3,6) from this traveling wave to 
the M = 2 level of the hierarchy, but does not predict the three M = 3 families of solutions 
that connect this traveling wave to interior bifurcations on paths B, C and D, nor the 
M = 4 solution connecting this traveling wave to path A from the previous example. In 
Figure |5j we show six solutions on this M = 4 path labeled E in the bifurcation diagram. 
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Figure 4: Degenerate bifurcation from M = 1 to M = 2,3,4 with ao held fixed. 

particle trajectories particle trajectories particle trajectories 

0.75 f~ ~~ If ' ' 1 : 




-0.5 0.5 -0.5 0.5 -0.5 



Figure 5: Particle trajectories along path E in the bifurcation diagram of Figure |4| 
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Figure 6: Left: Bifurcation diagram showing a path of time-periodic solutions at level 
M = 3 connecting two M = 2 solutions. Right: Three solutions on this path. 

When 9 = + , (where path E meets path A), a fourth particle nucleates at the origin 
without destroying periodicity of the other three. As 9 increases to 1, the trajectories of the 
particles j3j (t) undergo several topological changes (that determine which particles exchange 
positions over one period) until the trajectories coalesce into a single circular orbit at the 
degenerate traveling wave (with each particle moving counter-clockwise one and a half times 
per period). 

Finally, in Figure |6j we show a path of solutions at level M = 3 that connects two 
M = 2 solutions by interior (non-degenerate) bifurcations. The parameters of this path are 

AT = 3, £=(1,1,1), v = (-2,0,1). (78) 

We parametrize this path by 

*i = ^> 52 = S3 = ^T' (0<#<1)- (79) 

When 9 = 0, we obtain the M = 2 family of solutions that bifurcates from a two-hump, 
right traveling wave with indices (2, 1, 1, 1). For this family of solutions, the mean is related 
to the parameters s% = S2 and §2 = S3 via oq = 3 — 2s2', hence, holding the mean fixed 
requires that S2 remains constant. We will not reach the one-hump, right traveling wave at 



the bifurcation (1, 1, 1, 1) unless we increase oq to 3. As we increase 9 in (79), the trajectory 
of the particle that nucleates at the origin at 9 = + grows and merges with the trajectories 
of the original two particles through three topological changes: one at Fl, one not shown, 
and one at F3. Eventually, when 9 = 1, path F joins path G connecting (2,1,2,5) to 
(3,-1,1,5). Since the two-hump traveling wave moves to the right while the three-hump 
traveling wave moves to the left, the solution at F3 involves particles moving clockwise for 
part of their orbit and counter-clockwise at other times, leading to an interesting three- 
particle trajectory with 5-fold symmetry. 
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A Zeros of the polynomial P(-,X) 

As mentioned in the proof of Theorem [2j there is a small gap in the paper [22] showing 



that the multiperiodic solutions (46) satisfy the Benjamin-Ono equation, for the bilinear 



formalism used to derive these solutions requires that the zeros of / (or /') in (47) and 



(48) lie in the upper (or lower) half of the complex plane. In this appendix, we prove the 
equivalent assertion (in the case that all the k% are integers) that the zeros of the polynomial 



P(-,X) in (43) lie inside the unit disk. The two key ideas of this proof, namely showing 



that P (or /) has a representation as a determinant, and that the matrix is non-singular for 
\z\ > 1, are essentially due to Matsuno [16] and Dobrokhotov/Krichever |10j . respectively. 

Theorem 3 Suppose M > 2, k\, . . . , fejvf € N, v\, . . . , vm G R, x w , ...,xmo G % and 

< q\ < pi < g 2 < P2 < ■ ■ ■ < qM < VM- (80) 

Let M = {1, . . . , M} and define the polynomial 



p( Z , a) = Yi c s>? s 'z ks , c s =i n n bi , 

s&v(M) \(i,j)eSxS' / \jes' / 

where S' = M\S, k s = ^ i&s k-i, V& = Y^jeS' v v and 



(81) 



(Pj ~ ?»)(© "ft) 



bj 



Qj o -ikjXj 



j (<ij - ii)(pj - Pi) v Pj 

Then all the zeros Pi of P(-,X) lie inside the unit disk A G C provided \X\ = 1. 
PROOF: First we show that 



(82) 



P(z,X) 



M 



M 



-j=i J L i=i 

where the M x M matrix R(z, X) has entries 

1 

Rij(z, X) = n(z, X)5ij H — , 

Pi Qj 



(M) 



det R(z, A), 



(83) 



n(z, A) 



b7 iX-»izi« M 



Pi ~ Qi 



IK 1 - 



(84) 



j¥=i 



The symbol Yif^i indicates a product over j £ M omitting j = i, while n!<j^ i s a product 
over all pairs G M 2 such that i < j. By expanding det R = ^ CT sgn(<r)i2j )CT (j) and 

collecting like products of the r,(z, A), we find that 



det R = ^2 Y[r 

SeV(M) L ieS J 



deti?^', (Rs' 



1 



Ps' ~ Qs' 



(85) 
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Here {S[, . . . , S' n } is an enumeration of S' , i.e. Rs' is the nxn (Cauchy) matrix obtained by 
removing the rows and columns with indices in S from the Cauchy matrix {{pi — qj)" 1 }^^, 
and det R& is taken to be 1. The determinant of a Cauchy matrix is well-known, giving 



deti? 



S' 



Ui<j(Pi-Pj)((lj ~ ft) 



n 



, Pi ~ Qj 



(5') 



(86) 



Thus, the right hand side of (83) is equal to 



M 



■(£') 



i<j 



(M) 



e ii • ii mi ii ii 



seV(M) L jeS' 



-is 5 j^i 



i<j 



(87) 



If i and j are both in S or both in S' , the terms a?- in the final product cancel with 
corresponding terms in one of the previous two products. If i E S and j G 5", one of the 
factors of ajj (or a,-, if i > j) in the final product cancels with a^ 1 in the middle product, 
leaving behind Yl(i j)esxs' a *i' as required. 

Next we show that R(z, A) is invertible for \z\ > 1 and |A| = 1. Fix such a z and A. 
Define eL = \b7 l \~ Ui z kl \ so that 



M 

n 



(Pi - <?;) 2 (Pj - Qi)(qj ~ Pi) 



di > 1. 



Suppose for the sake of contradiction that there is a non-zero vector 7 £ C M such that 
i?7 = 0. This means that 



X] ( + ) 7j = HTi + ^(Pi) =0, (1 < i < M), i;{k) := ^ 

f=i V Pi q i J 7=1 



7j 



Then we define 



and observe that 

res k=Pi <p + res k=qi 



E 



7i 



E 



n 



7» 



P< 



n 



1j 



n ~ Pj 



(d 2 - 1) > 0, (1 < i < M) (90) 



where the inequality is strict if % 7^ and we used \ip(pi)\ 2 = | r i| 2 |7i| 2 - Thus, the sum of 
all the residues of cj>(k) is strictly positive, contradicting (p(k) = 0(k~ 2 ) as k — > 00. □ 
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B Uniqueness of Periodic Traveling Waves 



In this section we elaborate on the paper [5] showing that the only traveling wave solutions of 
the Benjamin-Ono equation are the solitary and periodic wave solutions found by Benjamin 
in [TJ. The purpose of this section is to modify their argument to prove that we have found 
all 27r-periodic traveling solutions, and to simplify part of their analysis. 

Consider any non-constant, 27r-periodic traveling solution of Q. After a transformation 
of the form (22), we may assume u{x) is a stationary solution satisfying 

uu x = Hu xx , (x G M/2vrZ). (91) 

After a translation, we may assume u x (0) = 0. Integrating once, there is a constant c such 
that 

1 c 2 

-u 2 = Hu x + -, (c>0). (92) 

The integration constant must be positive since Hu x is the derivative of a periodic function 
while the left hand side is positive. Now define the holomorphic function 

/i(*) = Z r^K&> (z £ A). (93) 



7T Jq K" — 6 

A direct calculation using Fourier series shows that 

f 1 (z + ) = a + u{9) + iHu{9), z = e ie , a = — [ u{9)d9, (94) 

where f\{z + ) is the limit of /i(C) as C approaches z from the inside. Now define 

hiz) = fi(e iz ) - ao, lmz>0. (95) 

Then f2{z) is analytic and bounded in the upper half-plane and satisfies 

f 2 (x) = u(x) + iHu(x), x€R. (96) 

Next we extend u{x) to the upper half-plane via u(x,y) = Re{/2(x + iy)} and define 

U(x,y) = u{x/C ^ + C , (97) 
2c 



where c was determined by u(x) in (92). We then have 

u y (x,0) =Re{i&(x)} = -Hu x {x) = ° ~™ , U y (x,0) = U(x,0) - U(x,0) 2 . 

Amick and Toland |5 showed that any non-constant, bounded, harmonic function U(x, y) 
defined in the upper half-plane and satisfying the nonlinear Neumann boundary condition 
U y = U — U 2 on the real axis as well as U x (0, 0) = 0, is given by 

U(x,y) = Re{f(x + iy)}, (y > 0), (98) 
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where f{z) is the (unique) solution of the complex ordinary differential equation 

/(0) = 17(0,0), a= y/2U(0,0)- C/(0,0) 2 (99) 



over the upper half-plane. Moreover, such a solution U(x,y) will satisfy {7(0, 0) S (0, 1) U 
(1,2], and the case £7(0, 0) = 2 corresponds to the solitary wave solution U(x, 0) = 2/(l+x 2 ), 
which is ruled out by the assumption that u in (91 ) is periodic. Rather than treat the cases 
£7(0,0) E (0,1) and {7(0,0) £ (1,2) separately as was done in [5], we choose the unique 
/3 G (-1, 1) such that 

(1 + (3) 2 1 - f3 2 



1/(0,0) 



1 + /3 2 



1 + /3 2 



and check directly that 



/(*) = o 



1 + 



2(3 







(100) 



:ioi: 



satisfies (99). Since U(x,y) = He{f(x + iy)} is x-periodic with (smallest) period while 
u(x,y) = [2cU(cx,cy) — c] is x-periodic with period 27r, it must be the case that c = N/a 
for some positive integer N. But then 



u(x) = Re{Nf 3 (Nx)}, h(z) 



2 /P 
a Va 



1 _ 1 - 3/? 2 4/3e iz 
a ~ 1 -P 2 + 1 — jSe* 



(102) 



i.e. ii(x) is one of the A^-hump stationary solutions discussed in Section [3] 
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